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ABSTRACT 

We present the first detection of a rotationally affected series consisting of 36 consecutive high-order 
sectoral dipole gravity modes in a slowly pulsating B (SPB) star. The results are based on the 
analysis of four years of virtually uninterrupted photometric data assembled with the Kepler Mission, 
and high-resolution spectra acquired using the HERMES spectrograph at the 1.2 meter Mercator 
Telescope. The spectroscopic measurements place KIC 7760680 inside the SPB instability strip, near 
the cool edge, given its fundamental parameters of T e g = 11650 ± 210 K, logg = 3.97 ± 0.08 dex, 
microturbulent velocity = O.O^q'q kms -1 , v sin* = 61.5 ± 5.0kms -1 , and [ M/H ] = 0.14 ± 0.09 dex. 

The photometric analysis reveals the longest unambiguous series of gravity modes of the same degree 
l with consecutive radial order n , that carries clear signatures of chemical mixing and rotation. With 
such exceptional observational constraints, this star should be considered as the Rosetta Stone of 
SPBs for future modelling, and bring us a step closer to the much needed seismic calibration of stellar 
structure models of massive stars. 

Keywords: asteroseismology — stars: variables: general — stars: individual (KIC 7760680) — stars: 
fundamental parameters — stars: oscillations — stars: rotation 


1. INTRODUCTION 

Slowly pulsating B stars (SPB stars) are non-radial 
multi-periodic oscillators on the main sequence between 
B3 and B9 in spectral type, 11 000 K and 22 000 K in ef¬ 
fective tempera ture, and 2.5 Mq and 8 Mq in mass (e.g., 
lAerts et~aLll2010l . Chapter 2). They pulsate in high-order 
gravity modes (with periods typically from 0.5 to 3 days) 
which are driven by the K-mechanism operating due to 
the op acity-bump associated with the iron-group ele - 
ments (lDziembowski et aLlll993tlGautschv fc Saiolll993 ). 
The g modes of the same degree i with consecutive ra¬ 
dial order n are expected to be almost equally spaced in 
period, and deviations from this equal spacing carry in- 
formati on a bo ut the physi cal processes in the near-core 
region dMiglio et all I2008t ). This allows us to use as¬ 
teroseismology to detect and interpret the fingerprints of 
mixing processes, such as core overshoot, diffusion, or ro¬ 
tationally induced mixing (originating from, e.g., internal 
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differential rotation or meridional circulation), that influ¬ 
ence the stellar lifetime significantly. Our understanding 
of these processes and their dependencies on different pa¬ 
rameters (such as mass, rotation rates, or the presence of 
a magnetic field) is highly unsatisfactory to provide the 
much needed calibration for stellar structure and evolu¬ 
tion models of the more massive stars with convective 
cores and radiative envelopes. This is partly due to the 
low number of in-depth studies of SPB stars, partly due 
to the low number (in terms of what is required for a 
secure forward modelling) o f identified modes in the ma¬ 
jority of these studies ('e.g.. !Aertsll2015il . Since the men¬ 
tioned models are cornerstones in many fields of stellar 
and galactic astrophysics, these shortcomings must be 
remedied. 

The Kepler M ission does not only excel in exoplanet 
detections (IBorucki et al.||2 010D. but it is al so a goldmine 
for asteroseismology (lGilliland et all 12010( 1. Thanks to 
the micromagnitude precision and the four year time- 
base of the initial mission, previously inaccessible diag¬ 
nostic methods are now within reach. This led to the 
fist actual seismic m odelling of an SPB star presented 
bv iPaoics et al. • : 201 li i. thanks to the unambiguous de¬ 
tection of 19 consecutive rotationally split i = 1 grav- 
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ity modes in the Kepler photometry of the slowly ro¬ 
tating star KIC 10526294. The detection of such pe¬ 
riod series combined with the approximate knowledge of 
fundamental parameters eliminates the time-consuming 
and resource-intensive need of individual mode identi¬ 
fication from follow-up multi-colo ur photometry and /or 
line-profile variation studies fe.g.. lDe Cat et aLll2005l i. 

KIC 7760680 {my = 10.3) was selected as a Kepler 
GO Program target together with KIC 10 526294 an d 6 
other SPBs (the procedure is described in [Panics et al.l 
I2013T ) . In Section[5] we show that the star is very close to 
KIC 10526294 in terms of the fundamental parameters, 
but it is rotating faster, which gives us the opportunity 
to study the effects of rotation. In Section[3] the results 
from the Kepler photometry are presented, specifically 
the detection of a rotationally influenced long period se¬ 
ries. In Section0]we discuss the implications and future 
prospects of this detection. 

2. SPECTROSCOPY 

To confirm its initial classification, KIC 7760680 was 
obser ved using the Hermes spectrograph dRaskin et al.1 
1201 If) installed on the 1.2 metre Mercator telescope on 
La Palma (Spain). A total of 8 spectra were taken: four 
consecutive exposures of 1000 —1800seconds in 2010, one 
900 s exposure in 2012, and three 1800 s exposures - each 
separated by two days - in 2014. The resolving power 
in the used HRF (high resolution fibre) mode was R ~ 
85 000, while the achieved signal-to-noise ratio (S/N) per 
exposure was typically ~ 30 in 2010 and 2012, and ~ 60 
in 2014. 

The raw exposures were reduced using the dedicated 
HERMES pipeline, and we worked with the extracted 
cosmic-removed, log-resampled (which includes barycen- 
tric correction), merged spectra afterwards. As a last 
reduction step, all s pectra were r e ctified following the 
process described bv iPapics et al.1 (|2013f ). We did not 
find any significant radial velocity variations in the data 
up to a level of ~ 5kms _1 that could indicate bina- 
rity, thus we immediately proceeded with calculating an 
average spectrum (using the inverse of the variance as 
weights). This mean spectrum with a S/N of ~ 120 was 
used to determine the fundamental para m eters. _ 

We used the GSSP co de package (jLehmann et al.l[2011b 
iTkachenko et al.1 [2012) to determine atmospheric pa¬ 
rameters of the star, such as the effective temperature 
(T e ff), surface gravity (log g), microturbulent and pro¬ 
jected rotational velocities (£t and z;sinz, respectively), 
and metallicity ([M/H]). The method is based on com¬ 
parison of the observations with a grid of synthetic spec¬ 
tra; a x 2 merit function is used to judge on the quality 
of the fit. The error bars are represented by 1 -a con¬ 
fidence levels that were computed from a y 2 statistics, 
and include all possible correlations between the param¬ 
eters and possible imperfections in continuum normal¬ 
ization, but exclude a ny other possible uncertainties such 
as NL TE effects (e.g JPuls et alJl2005l : iNieva fc Przvbillal 
2012) and/or incorrect atomic data. 

The spectrum was fitted in the wavelength range be¬ 
tween 4250 and 5800 A. The bluer part of the spectrum 
was omitted due to increasing uncertainties in the con¬ 
tinuum normalization and decreasing S/N of the data, 
whereas the red part was excluded due to contamina¬ 
tion by telluric lines. After the first run of the code, it 



Figure 1. Kiel diagram of a sample of B-type stars near the 
main sequence - for which an in-depth seismic analysis was carried 
out - observed by CoRoT (circles) and by Kepler (squares). The 
SPBs with an observed period series are plotted using empty sym¬ 
bols. The dot-dashed line represents the zero-age main sequence 
(ZAMS), while the dashed line the terminal-age main sequence 
(TAMS). The thin gray lines denote evolutionary tracks for selected 
masses, while the SPB (gray area) and /3 Cep (cool edge plotted us¬ 
ing a dashed line) instability strips are also shown. These were all 
calculated for Z = 0.015, X = 0. 7, using OP opacitie s (ISeatonj 
120051) . A04 heavy element mixture IIAsplund et al j|2005f ). and the 
linear nonadiabatic code by Dzicmbowski (1977), with the details 
of computations as in Pamyatnykh (1999). Different error bars 
reflect differences in data quality and methodological approach. 

already became clear that the fit is insensitive to the vari¬ 
ations of microturbulent velocity: our best fit suggested 
£t of 0.0 km s -1 with a 1-cr uncertainty level of 2.0 km s _1 . 
Moreover, we found the microturbulent velocity to corre¬ 
late strongly with the metallicity of the star, which varied 
between the solar value (0.0 dex) and 0.15 dex, for the 
range from 0.0 to 2.0kms _ . Fortunately, He lines are 
particularly sensitive to the variations of microturbulent 
velocity in B-type stars, and give little to no contribution 
to the overall stellar metallicity. Thus, we used helium 
lines to determine the microturbulent velocity parameter 
to be O.Olo'gkms - , and fixed it in the next step while 
fine-tuning the remaining four atmospheric parameters. 
This procedure was repeated several times until conver¬ 
gence was achieved; the 1-cr uncertainties in were prop¬ 
agated when computing the error bars in the other four 
parameters to take into account possible correlations. 

The derived parameters are T e s = 11650 ± 210 K, 
log g = 3.97 ± 0.08 dex, £ t = O.Olg;® kms -1 , v sin i = 

62±5kms _1 , and [M/H] = 0.14±0.09dex. They place 
the star inside the SPB instability strip (see Figure[l]) 
with a spectral type of B8V (based on T e e and logg 
values by using an i nterpolation in the tables given by 
iSchmidt-Kaleil 19821) . The obtained model fit to the 
mean spectrum is shown on Figure[2j We point out that 
the relatively large uncertainty of usinz reflects the ef¬ 
fect of un resolve d pulsational broadening due to g-modes 
(e.g.. lAerts et al.ll2014D . 

3. PHOTOMETRIC OBSERVATIONS AND FREQUENCY 
ANALYSIS 

KIC 7760680 was observed for a total of 1470.5 days 
(from the first commissioning quarter Q0 until the end of 
the last science quarter Q17, ~ 4 years) in Long Cadence 
(LC) mode with a cadence of 29.43 minutes by the Kepler 
satellite. We constructed light curves for each quarter 
from the target pixel Hies using custom masks. We added 
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Figure 2. Selected regions of the observed mean HERMES spec¬ 
trum (gray) and the best fit model spectrum (black solid line). Two 
Balmer lines, and two wavelength ranges that are rich in metal lines 
are shown. A few stronger lines are labelled, and an offset of 0.15 
units in flux are introduced between the selected regions for clarity. 

extra pixels with significant flux present beyond the ex¬ 
tent of the default mask, which resulted in quarterly 
curves with less instrumental trends tha n the standar d 
extraction (for more details see, e.g.. [Panics et al.ll2013l) . 
There are no other contaminating sources in the vicin¬ 
ity. After correcting for_the inc orre ctly reported times in 
Q0-Q14 (as noted in lPanics et all[2014! ). we cleaned all 
quarters from clear outliers, and detrended them using 
a division with a second order (except for Q0, where we 
used a linear) polynomial fit. Then the counts were con¬ 
verted to ppm and quarters were merged to a continuous 
light curve with a duty cycle of 91.6% (part of which is 
shown on Figured]). 

We extracted the Fourier-parameters ( Aj amplitudes, 
fj frequencies, and 9j phases) of the pulsation modes 
following a standard iterative prewhitening procedure 
and applied a significance criterion to select the strictly 
significant o n es (se e Figured]), in the same way as in 
iPfipics et all (|2014l ). 

3.1. Combination Peaks 

The Scargle periodogram of KIC 7760680 shown in Fig¬ 
ured] shows clear groups of peaks. Following the ap¬ 
proach described bv lPanicsl ( 2012 1 we find that all peaks 
below 0.2 d -1 (2.31 pHz) and above 1.2 d -1 (13.9 pHz) 


are low order combination frequencies originating from 
independent parent modes in the range in between. Such 
harmonics and linear combinations are observed i n many 
classes of classical pulsators (see, e.g..lDolez et, al.ll2006l : 
iBreger et, all 1 20 111 : I Van R.eeth et, al.ll2015bH . and various 
non-linear mechan isms may b e responsible for their ap- 
peareance (IBreger et al.l 1201 lh . The spectrum of inde¬ 
pendent modes is simple, and shows a clear structure. 
No pressure modes are observed towards shorter periods, 
as expected from excitation calculations for this temper¬ 
ature range. 

3.2. A Series of Gravity Modes 

The frequency peaks in between the mentioned groups 
of combination frequencies fall in the expected range for 
high-order gravity modes of SPB stars. Here, we detect 
a clear series of 36 frequencies (listed in Table]]]) which 
follow a characteristic period spacing pattern (see Fig- 
ure|4]). There is a distinct cutoff in power towards longer 
periods at the frequency 0.63051(2) d -1 (7.2976(2) pHz), 
which we interpret as a good proxy of t he buoyancy cut - 
off frequency (e.g., iHansen et, alJll985t lTownsendll200f)f) . 
This cutoff value is higher than the one of KIC 10526294 
by 0.16d -1 (1.85/iHz), which must reflect a difference 
in mode energy leakage from the stellar atmosphere due 
to the difference in evolutionary state and/or rotation. 

I Townsend ( 2000 ) indeed concluded that the Coriolis force 
increases the cutoff frequency, although not so much for 
prograde sectoral modes. 

This is only the fourth main sequence B-type 
star where we find such a series of period spac- 
ing s (after HD 50230, HD 43317, and KIC 10526 294, 
by iDegroote et, al.l 120101 : iPaoics et, al.l 120121 12014 re¬ 
spectively). Two of them are ultra-slow rotators, 
KIC 7760680 is a moderate rotator, and HD 43317 rotates 
at 50% of its critical velocity. This new detection con¬ 
firms that gravity-mode period spacings occur in pulsat¬ 
ing B stars irrespective of their rotation rate and offers a 
great opportunity to probe the interior rotation pro file of 
such stars, as first tried for KIC 10526294 dTriana et all 

nil). 

The number of consecutive peaks in the series for 
KIC 7760680 is much higher than the previous record 
holder (19 peaks in KIC 10526294), and this is the first 
time that we detect a significant tilt and a long periodic 
pattern in the period spacing series. There is also an 
apparent bifurcation in mode amplitudes around 1.2 d, 
with every second peak being significantly stronger for a 
length of approximately 5 pairs. This implies that modes 
with even radial order get excited with different ampli¬ 
tude than modes with odd radial order near the maxi¬ 
mum power. As Figure|4] shows, this difference in am¬ 
plitude is not connected to mode trapping in the inner 
near-core regions; it rather seems to indicate a difference 
in the pulsation mode energy distribution in the envelope 
for odd versus even radial order dipole modes. 

4. DISCUSSION 

We have identified the clearest and longest period se¬ 
ries of consecutive high-order gravity modes in an SPB 
star to date. We see clear effects of rotation in the 
tilt of the period spacing series. Similar effects have 
also been found in y Do radus stars (|Bedding et al.l]20l4 
iVan Reeth et al.l l2015al Tbh. The series here is not only 
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Figure 3. One tenth of the reduced Kepler light curve (top panel, brighter points towards the top), and the distribution of significant 
frequencies (bottom panel, black vertical lines) in the frequency spectrum of KIC 7760680. Combination peaks are marked with gray dashed 
lines (darker color corresponds to simpler combinations). The region of maximum power is shown in detail on Figure[4] 


tilted, but also shows periodic oscillations - with an 
unprecedented clarity - which carry information about 
the size of the convective core (the age of the star) 
and the mixing processes inside the star. Unlike for 
KIC 10526294, which turned out to be too young to show 
these oscillations, or the first two detections where the 
period spacing series were too short to be conclusive, fu¬ 
ture detailed seismic modelling of this star will be able to 
put stringent constraints on the age, the core-overshoot 
parameter and diffusive chemical mixing properties of the 
star. 

A compatibility-check based on the expected period 
spacing behaviour relying on Ledoux splitting for high- 
order g-modes (see, e.g., Chapter 3 in lAerts et al.ll2010h 
is shown in the bottom panel of Figure[U The overall 
trend in the observed period series (both the tilt and the 
periodic deviations from it) can be reproduced from an 
appropriate stellar model with M. = 3.3 Af®, X = 0.71, 
Z = 0.02, central hydrogen fraction X c = 0.44, and 
core overshoot a ov = 0.3, using t = l,m = 1 inodes 
assuming a rotation rate that is compatible with the 
measured wsini. The rotation rate might place the ob¬ 
served modes in the gravito-inertial regime, where the 
per turb ative treatm ent of rotation is not valid anymore 
ijBallot et al.l 1 201 Of ) . Therefore, future in-depth mod¬ 
elling will have to take this into account after derivation 
of U(r), and also explain why the retrograde sectoral and 
the z onal d ipole inodes are not visible. It is noteworthy 
that iTownsendl (|2003D already provided quantitative ar¬ 
guments that favor the visibility of prograde dipole sec¬ 
toral modes, completely in agreement with our observa¬ 


tional results. Once we have a good seismic model of 
the star, we also plan to investigate the physical nature 
of the linear combination frequencies in terms of non¬ 
linear response of the stel l ar flux to one or a few high- 
amplitude inodes (iGarrido fe RodriguezllT 996i) or due to 
resonant mode coupling ( Buchler et al.lll997l) . The latter 
predicts specific resonances to occur among the various 
heat-driven inodes. In line with such theoretical predic¬ 
tions, we can then check if there is more a posteriori 
information in the frequency spectrum beyond the de¬ 
tected i = 1 series. Such further interpretation can only 
be done after the dominant dipole series is adequately 
modelled. 

This star has the potential of becoming the seismic 
Rosetta Stone of SPB stars. While ground based mode 
identification efforts can at best only deli ver constraints 
on a few osc illation freq uencies (e . g.. iChanellier et al.1 
l2000t iDaszvriska-Daszkiewicz et al.l 120151) . the unam¬ 
biguously detected series of 36 modes carry the seismic 
information needed to understand the physical processes 
at play inside SPB stars better. For the near future, we 
need to find more SPBs with such series, and we plan to 
invest in their in-depth seismic modelling. 

The research leading to these results was based on 
funding from the Fund for Scientific Research of Flanders 
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Council of KU Leuven under grant GOA/2013/012, Bel¬ 
gium. This research was also supported in part by the 
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1.1 1.2 
Period (d) 


Figure 4. Zoom-in of the period spectrum of KIC 7760680 (top panel) with the members of the period spacing series marked using vertical 
dashed lines. The observed (middle panel) period spacing is shown along with a series from an appropriate stellar model, discussed in the 
text, as comparison (bottom panel). 


National Science Foundation of the United States under 
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Table 1 

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the detected period series. The displayed S/N values 
are calculated in a window of 1 <1 “ 1 in the residual spectrum centered on the given frequency. Numbers in parentheses are the formal 

errors of the last significant digit. 


# 

P 

d 

/ 

d" 1 

J 

ppm 

0 

27r/rad 

S/N 

1 

0.86930(2) 

1.15035(3) 

0(7) 

-0.43(7) 

4.2 

2 

0.90147(1) 

1.10930(1) 

40(2) 

-0.21(4) 

13.3 

3 

0.93036(1) 

1.07485(1) 

44(2) 

0.49(4) 

8.5 

4 

0.95674(1) 

1.04521(1) 

45(2) 

0.42(4) 

8.6 

5 

0.981691(6) 

1.018650(6) 

219(3) 

0.06(2) 

22.0 

6 

1.007872(2) 

0.992189(2) 

1386(7) 

-0.074(5) 

60.3 

7 

1.033219(4) 

0.967849(4) 

526(5) 

0.28(1) 

35.1 

8 

1.05421(1) 

0.94858(1) 

78(2) 

0.48(3) 

11.4 

9 

1.073164(5) 

0.931824(4) 

620(6) 

0.215(9) 

36.6 

10 

1.094833(4) 

0.913381(3) 

921(6) 

0.399(7) 

48.6 

11 

1.117507(1) 

0.894849(1) 

3400(11) 

0.326(3) 

84.8 

12 

1.139980(8) 

0.877208(6) 

225(3) 

0.33(1) 

21.7 

13 

1.158919(1) 

0.862873(1) 

9876(26) 

0.409(3) 

96.9 

14 

1.175637(7) 

0.850603(5) 

255(4) 

-0.37(1) 

24.2 

15 

1.192566(1) 

0.838528(1) 

5017(18) 

-0.201(4) 

86.4 

16 

1.210189(7) 

0.826317(5) 

256(4) 

0.25(1) 

23.4 

17 

1.228020(3) 

0.814319(2) 

1755(9) 

-0.298(5) 

62.4 

18 

1.24544(1) 

0.802928(9) 

130(3) 

0.36(2) 

13.4 

19 

1.259035(6) 

0.794259(4) 

464(5) 

0.14(1) 

32.3 

20 

1.27561(2) 

0.78394(1) 

60(2) 

0.44(3) 

9.8 

21 

1.28944(3) 

0.77553(2) 

31(1) 

-0.37(4) 

6.7 

22 

1.30351(3) 

0.76716(2) 

39(2) 

-0.16(4) 

7.8 

23 

1.31591(2) 

0.75993(1) 

49(2) 

0.34(4) 

8.0 

24 

1.32771(1) 
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